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We have directly imaged the anisotropic nonlinear Meissner effect in an unconventional super- 
conductor through the nonlinear electrodynamic response of both (bulk) gap nodes and (surface) 
Andreev bound states. A superconducting thin film is patterned into a compact self-resonant spiral 
structure, excited near resonance in the radio-frequency range, and scanned with a focused laser 
beam perturbation. At low temperatures, direction-dependent nonlinearities in the reactive and 
resistive properties of the resonator create photoresponse that maps out the directions of nodes, 
or of bound states associated with these nodes, on the Fermi surface of the superconductor. The 
method is demonstrated on the nodal superconductor YBa2Cu307_a and the results are consistent 
with theoretical predictions for the bulk and surface contributions. 



Introduction - The Meissner effect is the spontaneous 
exclusion of magnetic flux from the bulk of a supercon- 
ductor. In the presence of a magnetic field, a supercon- 
ductor must invest kinetic energy in a supercurrent flow 
to screen out the applied field. This reduces the free en- 
ergy difference between the superconducting and normal 
states, resulting in a reduction in magnitude of the su- 
perconducting order parameter. This in turn changes the 
properties of the superconductor (e.g. magnetic penetra- 
tion depth, diamagnetic moment, etc.) and is referred 
to as the nonlinear Meissner effect (NLME). Conven- 
tional (fully gapped) superconductors show the strongest 
nonlinearities near T c , and have exponentially suppressed 
nonlinear response at low temperatures, T -C T c . Uncon- 
ventional superconductors with nodes in the supercon- 
ducting energy gap are expected to have a strong nonlin- 
ear Meissner effect at low temperatures, due to the nodal 
excitations out of the superconducting ground state [lj. 
In addition this nonlinear response should be anisotropic, 
reflecting the locations of nodes of the gap on the Fermi 
surface. 

For a d X 2- y i gap on a circular cut of the cylindrical 
Fermi surface, typical of cuprate superconductors, the- 
oretical treatments of the anisotropic NLME (aNLME) 
predicted a 1 / y/2 anisotropy at zero temperature [2, E| . 
Later, the theory was re-formulated in terms of nonlin- 
ear microwave intermodulation response of a nodal su- 
perconductor and the temperature dependence of 
the NLME and its anisotropy was worked out. The de- 
pendence of the superfluid density n s on supercurrent 
density J s , for sufficiently small currents, is given by 
n s (T,J s ) = n s (T)[l - b e (T)(J s / J c )% where J c is the 
critical current density, and &e (T) is the nonlinearity co- 
efficient dependent on the direction (angle 9) of J s 0. It 
was found that the anisotropy in the NLME of cuprates 



is weak at high temperatures, and only becomes signif- 
icant for T/T c < 0.6 0. In addition, it was found that 
b e (T) is expected to grow as l/T for T/T c < 0.2 0, 
before crossing over to another temperature dependence, 
depending on the purity of the material 0, [H Q . 

Early experiments to detect the aNLME of cuprates 
through transverse magnetization 0, Hj], magnetic pen- 
etration depth , and surface impedance [12j mea- 
surements did not establish conclusive evidence of the 
effect Later, sensitive nonlinear microwave measure- 
ment techniques established the existence of the NLME 
in cuprates from the temperature dependence of the in- 
termodulation power at low temperatures 13 -llj], al- 
though the anisotropy of the effect was not measured. Al- 
though other thermodynamic measurements have shown 
evidence of an anisotropic energy gap in various super- 
conductors 

EMI, there are no unambiguous measure- 
ments of the aNLME in cuprates, to our knowledge. In 
contrast to quasiparticle transport methods Til- 18j. the 



present method utilizes the bulk superfluid response to 
identify the nodal directions [l~9j |. 

An additional contribution to the NLME in uncon- 
ventional superconductors arises from Andreev bound 
states (ABS) [20j, created, for example, on the (110) 
surface of a d X 2- y 2 superconductor. These states giv e 
rise to a paramagnetic contribution to the screening [21| . 
For cuprates, (110) interfaces occur at twin boundaries, 
which are formed spontaneously during epitaxial film 
growth. The NLME associated with ABS has been estab- 
lished by tunneling [22j , and penetration depth measure- 
ments @, [24| , for example. Theory by Zare, Dahm and 
Schopohl (ZDS) predicts an aNLME associated with ABS 
having a strong temperature dependence at low temper- 
atures, eventually dominating that due to nodal excita- 
tions dn. 
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Here we establish a new method to both quantitatively 
measure and image the aNLME from nodes in the su- 
perconducting gap using a novel RF resonant technique 
combined with laser scanning microscopy. 

Experiment - We employ a self-resonant superconduct- 
ing structure based on a thin him Archimedean spiral 
geometry. The spiral has an inner diameter of 4.4 mm, 
an outer diameter of 6 mm, and consists of 40 turns of 
nominally 10 /im width thin him stripe with 10 /im spac- 
ing, winding continuously from the inner to outer radii 
with Archimedean form (see the schematic in Fig. 1(a)). 
The structure has considerable inductance and capaci- 
tance, making it a compact self-resonant meta-atom for 



use in superconducting metamaterials [261 . I27j . When 
made from superconducting materials such as Nb and 
YBa2Cu307_,5 (YBCO), the spiral has a fundamental 
resonance in the vicinity of 75 MHz, and many overtones 
extending above 1 GHz. The distribution of standing 
wave currents on the spiral in the hrst few modes are 
well approximated as those of a resonant vibrating string 
held hxed at both ends, and then wrapped into a spiral, 
as verified b y d etailed laser scanning microscope (LSM) 



imaging 
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A unique property of the resonant spi- 
ral is the fact that currents flow in a uniform manner 
around the Fermi surface of the superconductor. 

The epitaxial and in-plane oriented YBCO films on 
LaA103 (LAO) were deposited by in-situ off-axis mag- 



netron sputtering to a thickness of 300 nm 31] . Similar 
YBCO films on Ce02 buffered sapphire and on MgO were 
deposited by thermal coevaporation to a thickness of 300 
nm and 700 nm, respectively [32j]- The 200 nm thick 
Nb films were deposited on quartz substrates at room 
temperature by RF sputtering [2|| 28 1. All films were 



patterned into spiral resonators by contact photolithog- 
raphy and either wet or dry chemical etching. 

The LSM operates by scanning a focused laser spot 
(wavelength 640 nm) over the surface of the spiral while 
it is excited near its RF resonance with excitation and 
pickup loops above and below the plane of the film. De- 
tails of the RF and laser excitations of the spirals have 
been previously published 26, 28- 3(|. The spiral devel- 
ops a photoresponse due to pair-breaking and localized 
heating, producing a change in its resonant frequency 
and quality factor H, 34|. The laser intensity is mod- 
ulated at 100 kHz, and the changes in RF transmission 
at a fixed frequency are phase sensitively detected, re- 
sulting in a photoresponse (PR) signal consisting of both 
a magnitude and phase. The peak laser intensity can 
be varied between 150 nW and 1.6 mW. The laser spot 
(about 20 /im diameter) is scanned over the sample to 
create a two-dimensional PR image. 

The aNLME is expected to create the following PR 
contrast [H, EIJ. For a d X 2- V 2 gap on a circular Fermi 
surface one expects PR ~ 1 + sin 2 (26) with 9 = along 
the Cu-O-Cu bond direction of YBCO (see Fig. 1(a)). 
In addition, the PR amplitude should vary as 1/T 2 for 
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FIG. 1. a) Schematic diagram of spiral geometry, definition 
of radial (p) and angular (0) coordinates, and directions of 
the crystallographic a- and b-directions, along with the orien- 
tation of the d-wave gap in YBCO. b)-d) LSM photoresponse 
images of a YBCO/LAO spiral resonator in the third har- 
monic mode at temperatures of b) 81.9 K, c) 4.40 K, d) 3.90 
K. Also shown is the d-wave gap orientation as determined 
from twin domain boundary directions in the substrate. 



T/T c < 0.2 and J RF /J C < T/T c , where J RF is the RF 
current density induced in the sample. It is also expected 
that the LSM PR is proportional to RF power and pro- 
portional to the intensity of the laser light. As deduced 
from ZDS, the contribution to LSM PR from ABS should 
be centered on (110) surfaces of the material, have a 
stronger temperature dependent PR ~ 1/T 4 , have the 
opposite sign of PR from the nodal aNLME (due to the 
paramagnetic nature of the ABS), and be a linear func- 
tion of RF and laser power up to RF fields on the order 
of a few mT [H. 

Results - Figure 1 shows LSM PR images of a single 
YBCO/LAO spiral in a resonant mode at three different 
temperatures, 81.9 K, 4.4 K and 3.9 K. The image is of 
the third harmonic mode and shows three nearly circu- 
lar distinct bands of enhanced PR, corresponding to the 
three half- wavelengths of current distributed over the 40 
turns of superconducting wire 



28-30!. One notes that 



the PR is isotropic in its angular distribution around the 
spiral at 81.9 K (Fig. 1(b)). However, the LSM PR image 
of the same mode at 4.40 K (Fig. 1(c)) shows four distinct 
enhancements of PR at regular angular intervals around 
the spiral. As the sample is cooled further (below about 
4.30 K for this particular sample), a new set of enhanced- 
PR features are observed in Fig. 1(d), but now rotated 
45° relative to those at 4.40 K. Images of Nb spirals of 
nominally identical geometry in the same mode do not 
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FIG. 2. Plot of radially-integrated and unwrapped photore- 
sponse (PR) vs. angle on the YBCO/LAO spiral at 81.9 K 
and 4.66 K, along with a fit of the low-temperature integrated 
PR(Q) to the simple d-wave angular dependence. Inset shows 
PR image of YBCO/LaA103 thin film spiral resonator taken 
at a temperature of 4.656 K in the 7 th mode at 450.22 MHz. 



show this anisotropic PR at any temperature below the 
transition temperature. The above anisotropy and rota- 
tion behavior is observed to hold for the YBCO/LAO spi- 
ral excited into all modes, from the fundamental through 
the 8*^ mode (the highest examined). 

Figure 2 shows the radially(p)-integrated PR as a func- 
tion of angle for a YBCO/LAO thin film spiral resonator 
taken at a temperature of 4.66 K in the 7 th mode at 
450.22 MHz, shown in the inset. The zero-angle (Cu-O- 
Cu bond) direction was determined from the visible twin- 
ning structure of the LAO substrate. Utilizing the fact 
that the YBCO film was grown epitaxially and in-plane 
oriented on LAO, the observation that the PR maxima 
are aligned along the linear twin domain structures in 
the substrate, and the fact that the twin boundaries are 
aligned at 45° to the Cu-O-Cu bond direction [36|, we 
deduced that the 8 = direction is at 45° to the twin 
domain structure of the substrate. Figure 2 shows a fit 
of the data to PR(<d) = 1 + J 4sin 2 (26), resulting in a fit 
value A = 1.05 ± 0.02, close to the value of 1 expected 
for a simple d X 2- V 2 gap on a circular Fermi surface. 

The PR images in the 'rotated' state at lower tem- 
peratures (Fig. 1(d)) also show 4-fold character, but 
are relatively diffuse in appearance, and display dramat- 
ically stronger PR as the temperature is lowered. To 
understand these properties, we examine the tempera- 
ture dependence and phase of the photoresponse along 
the nodal (6 = 7r/4) and anti-nodal (O = 0) directions 
of the superconducting gap. Figure 3 shows a plot of PR 
magnitude vs. temperature for nodal and anti-nodal re- 
gions of a YBCO/sapphire spiral resonator. In the nodal 
direction the PR is of significant magnitude at a tem- 
perature of 6.6 K, but decreases in magnitude as tem- 
perature is lowered. The PR goes to near-zero levels 
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FIG. 3. Plot of temperature dependence of photoresponse 
(PR) magnitude taken at nodal and anti-nodal locations in a 
YBCO /sapphire spiral resonator. The anti-nodal PR has op- 
posite sign to the nodal PR above the crossover temperature 
(about 5.5 K). Inset shows temperature dependence of nodal 
PR magnitude near its crossover temperature, along with a 
fit to the form predicted by ZDS. Other insets show LSM PR 
magnitude images of YBCO/sapphire below (4.405 K) and 
above (6.641 K) its crossover temperature at 546.555 MHz. 
The phase of the PR in these images differ by 7r, correspond- 
ing to a sign reversal of the response. 

and changes phase by ir radians at a point that we call 
the 'crossover temperature.' Below this temperature the 
nodal PR magnitude increases dramatically. In the anti- 
nodal directions the PR is found to be smaller than the 
nodal PR at 6.6 K, but increases monotonically in mag- 
nitude with decreasing temperature, and has the same 
phase as the nodal PR below the crossover temperature. 

A fit of the PR to the expected temperature depen- 
dence of nodal PR (temperature derivative of Eq. (16) 
of Ref. 25), namely PR ~ A/T 2 - B/T 4 (first term due 
to nodal NLME and second due to ABS NLME) [25j| . 
is shown in the inset of Fig. 3 for a YBCO/sapphire 
sample. The data is consistent with the predicted tem- 
perature dependence around the crossover temperature. 

In all of these experiments we have verified that the 
LSM PR in the nodal direction is proportional to both 
RF power and to laser intensity. The LSM PR in 
the anti-nodal direction is linear at low laser powers, 
but becomes strongly nonlinear at higher laser powers. 
YBCO films grown on other substrates (MgO and sap- 
phire) show the same general behavior described above 
as the samples on LAO. The insets of Fig. 3 shows im- 
ages of the YBCO/sapphire sample below (4.405 K) and 
above (6.641 K) the crossover temperature for this sam- 
ple. We find that the PR maxima are aligned along the 
nodal directions of YBCO at higher temperatures, and 
then show a 'rotation' by 45° at lower temperatures, al- 
though the crossover temperature is sample dependent 
(approximately 4 - 5 K for YBCO/LAO and 5 - 7 K 
for YBCO/MgO and YBCO/sapphire). This sample de- 
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pendence may be associated with different abundances of 
(110) surfaces and twinning in each sample. 

Discussion - A reorientation of specific heat oscilla- 
tions with magnetic field direction (qualitatively similar 
to the crossover observed here) has been reported in su- 
perconductors with anisotropic order parameters [13, HH . 
However, specific heat measures the low-energy quasi- 
particle density of states while our measurements 
probe the superfluid response and involve no DC mag- 
netic field or vortices. Therefore we do not believe that 
quasiparticle transport phenomena are the origin of the 
crossover observed in this paper. On the other hand, the 
observed sign difference of LSM PR between the nodal 
and anti-nodal directions above the crossover tempera- 
ture is consistent with the ZDS prediction for the NLME 
due to ABS [lH . The strong and monotonic temperature 
dependence of PR observed along the (110) directions 
of the film are also consistent with the ZDS prediction. 
The observed crossover temperature (4 K to 7 K) is on 
the order of that predicted by ZDS, namely Tc/k 1 ^ 2 ~ 
10 K, where k is the Ginzburg-Landau parameter. To- 
gether, these results confirm the basic predictions of ZDS 
for the relative contributions of the nodal and ABS con- 
tributions to the NLME. Note that the use of thin film 
superconducting spiral resonators that expose many dif- 
ferent crystallographic directions, and which have a large 
surface-to-volume ratio, accentuates the contribution of 
ABS to the photoresponse. 

In order to verify the above interpretation of the data, 
a number of potential artifacts and other interpretations 
of the PR images have been systematically explored. 
First is the possibility of anisotropic PR being created by 
the irregularly-shaped excitation and RF pickup loops. 
The anisotropic PR did not change when the excitation 
loop was rotated relative to the sample. In addition, 
the anisotropic PR pattern rotated when the sample was 
rotated relative to the loops. We do not observe a sig- 
nificant degree of PR anisotropy in Nb spirals measured 
under similar conditions. 

A second concern is that the PR temperature depen- 
dence arises from a growing thermal boundary resistance, 
Rbdy ~ 1/T 3 , between the thin film superconductor and 
the dielectric substrate [4(J • However the strong temper- 
ature dependence is not observed in Nb/quartz samples, 
eliminating most of the cryostat components from gener- 
ating this response. All YBCO samples (on LAO, MgO 
and sapphire) show the strong PR temperature depen- 
dence at low temperatures. Nevertheless, we do believe 
that a thermal blockade may set in for the YBCO/LAO 
samples at low temperatures (< 3 K), resulting in diffuse 
PR images. This interpretation was independently veri- 
fied through measurements of the PR thermal relaxation 
time as a function of temperature. 

A third concern is that the twinned structure of the 
LAO substrate traps heat and channels it in specific di- 
rections, giving rise to the observed anisotropy. How- 



ever, we observe the same anisotropy on YBCO/MgO 
and YBCO/sapphire films that do not have twinned sub- 
strates. The apparent rotation of the PR by 45° in a 
narrow temperature range in all of the YBCO samples 
also argues against a substrate-induced anisotropy. Next 
is the concern that the shape of the substrate dictates 
the anisotropy of the PR in the spiral. Numerical simula- 
tions of localized heating of a spiral on a square substrate 
under conditions similar to those of our YBCO samples 
showed a temperature anisotropy of only 1 part in 1000 
imposed by the shape of the substrate, whereas the ob- 
served anisotropy is on the order of 50%. Finally is the 
concern that defects in the lithographic pattern of the 
spiral (namely shorts between wires or opens) will create 
the anisotropic PR. We find that all YBCO spirals, from 
those with near-perfect lithography to those with defec- 
tive patterns, all show the anisotropy and temperature 
dependent properties described above. 

Conclusions - We have measured and imaged the 
anisotropic nonlinear Meissner effect in an unconven- 
tional superconductor arising from both bulk and sur- 
face mechanisms. Our results are in detailed agreement 
with aNLME theoretical predictions based on both nodal 
and ABS contributions to the photoresponse. This es- 
tablishes a new gap spectroscopy tool that is sensitive 
to nodes in the superconducting gap accessed by cur- 
rents flowing in the plane of the material, and can detect 
gap nodes directly, or indirectly through Andreev bound 
states, through their influence on the in-plane electrody- 
namics in the superconducting state. 
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